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Compact binary systems with neutron stars or black holes are one of the most promising sources 
for ground-based gravitational-wave detectors. Gravitational radiation encodes rich information 
about source physics; thus parameter estimation and model selection are crucial analysis steps for 
any detection candidate events. Detailed models of the anticipated waveforms enable inference on 
several jjaranieters, such as component masses, spins, sky location and distance, that are essential for 
new astrophysical studies of these sources. liowcvcr, accurate measurements of these parameters and 
discrimination of models describing the underlying physics are complicated by artifacts in the data, 
uncertainties in the waveform models and in the calibration of the detectors. Here we report such 
measurements on a selection of simulated signals added either in hardware or software to the data 
collected by the two LIGO instruments and the Virgo detector during their most recent joint science 
run, including a "blind injection" where the signal was not initially revealed to the collaboration. 
We exemplify the ability to extract information about the source physics on signals that cover the 
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neutron-star and black-hole binary parameter space over the component mass range 1 M0 — 25 Mq 
and the full range of spin parameters. The cases reported in this study provide a snap-shot of the 
status of parameter estimation in preparation for the operation of advanced detectors. 

PACS numbers: 04.80.Nn, 04.25. dg, 95.85.Sz, 97.80.-d 



I. INTRODUCTION 

General relativity predicts that binary systems of com- 
pact objects lose energy through the emission of gravita- 
tional radiation, a prediction confirmed through binary- 
pulsar observations [l]-(3]. During this process they emit 
a characteristic 'chirping' gravitational wave (GW) sig- 
nal of predominantly increasing amplitude and frequency. 
For neutron stars and black holes the signal enters the 
observational band of the initial (at 40 Hz) and advanced 
(at 10-20 Hz) ground-based laser-interferonietric detec- 
tors (e.g., [IHl]), sweeping through the detection band 
for a few seconds (depending on the masses of the ob- 
jects) until coalescence. 

The search for GW signatures of compact binary co- 
alescence in LIGO's most recent sixth science run and 
Virgo's science runs 2 & 3 is described in [7|[B]. No detec- 
tion was reported. However, the operation of advanced 
instruments - Advanced LIGO and Advanced Virgo - 
from 2015-1- [3 ^ [TO] suggests that the coalescence of 
compact binaries could be observed in the not too dis- 
tant future [11]. Once a detection candidate has been 
identified, the next step in the analysis is to measure the 
source parameters and test models that describe the un- 
derlying physics. This step is critical for enabling studies 
in astrophysics and fundamental physics in which GW 
observations are expected to provide a new view of rel- 
ativistic phenomena. One of the signals studied in this 
work represents an end-to-end test of this process: a de- 
tection candidate was identified with a false alarm rate of 
1 in 7000 years (1.4 x 10"'' yr""') [7] and was fully char- 
acterized in terms of its physical properties before it was 
revealed to be a "blind" hardware injection, i.e., a sim- 
ulated signal added by coherently actuating the mirrors 
of the LIGO and Virgo detectors, without the knowledge 
of the data analysts [T^ . 

The ability to accurately estimate the parameters of 
coalescing binaries, including the masses of the compo- 
nents, their spins, the location of a binary on the sky and 
its distance, and the implications for new insights into the 
underlying physical processes have been at the centre of 
several studies. Accurate measurements of the masses 
provide information about the mass distribution of (bi- 
nary) black holes and neutron stars, clues to determining 
the maximum mass of neutron stars [e.g.,[TH], the under- 
lying neutron star equation of state, the minimum mass 
of stellar-mass black holes and the presence or absence 
of the so-called "mass gap" [TH - fTO] . Spin measurements 
are a direct window onto the critical stage of common 
envelope evolution, and the details of the supernova pro- 
cesses and kicks. Localizing the merger of compact bi- 
nary systems on the sky may permit the identification of 



the host, and therefore an opportunity for in-depth stud- 
ies of the environments in which compact binaries form 
[e.g.. fT7H2T] . If the localization region in the sky is suf- 
ficiently small, a possible electro-magnetic counterpart 
could be found and multi-messenger studies of coalesc- 
ing compact objects would become possible [e.g.. l^2H^ . 
For instance, the question of whether compact binaries 
are the progenitors of short gamma-ray bursts could be 
definitively answered. The direct determination of the 
luminosity distance to the source opens new possibilities 
for low rcdshift cosmography. In general, coalescing bi- 
naries are a laboratory for tests of the behavior of gravity 
in the strong-field regime. Eventually, when multiple de- 
tections are available, studies of source populations will 
be possible [e.g., [551 as well as more stringent tests 
of the underlying source dynamics [IT] [55] . 

Most of these studies use theoretical estimates 
of parameter uncertainty based on the Cramer-Rao 
bound [29) . which should be valid in the limit of high 
signal-to-noise ratio (SNR). Initial detections may be too 
weak for this bound to provide useful guidance. There- 
fore, a complete Bayesian analysis like the one described 
below must be used to quantify parameter uncertain- 
ties. Other studies have relied on injections into syn- 
thetic data. In this paper, we will use injections into real 
data, which introduces a new set of challenges, including 
non-Gaussianity and non-stationarity. 

The fact that gravitational waveforms used in the anal- 
ysis are an approximation to the actual radiation pro- 
duced by astrophysical sources and that the measured 
strain is affected by the uncertainties in the instrument 
calibration [5Dl - [5^ represent additional challenges for 
making robust inference on the underlying physics. To 
study parameter estimation in this regime, we have ana- 
lyzed several artificial compact binary coalescence (CBC) 
signals added to real detector data, including the "blind" 
injection described above, added both in hardware and 
software to the data collected by the two LIGO instru- 
ments (Hanford and Livingston) and the Virgo detec- 
tor during the most recent joint science run, S6/VSR2-3. 
The use of injections has been, and continues to be, an 
essential means to validate the detection process, and 
as we report here, has been naturally extended to the 
source characterization stage of the analysis. Here we 
exemplify the ability to extract information about the 
source physics on a selected number of injections that 
cover the neutron star and black hole parameter space 
over the component mass range 1 Mq — 25 Mq and the 
full range of spin parameters. We consider a spectrum of 
realistic signal strengths, from candidates observed close 
to the detection threshold to high SNR events, and vari- 
ous relative strengths across the instruments of the net- 
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work. We analyzed the signals using a range of waveform 
models that demonstrate the interplay between (some) 
systematic bias and statistical uncertainty. To help val- 
idate our results, we carry out the analysis with several 
independent techniques; these are implemented within 
a specially-developed software package part of the LSC 
Algorithm Library, LALInference 33] . 

The paper is organized as follows. In Section |ll] we 
give a brief overview of the analysis method. While no 
detections were claimed in '7], simulated signals ('injec- 
tions') were added to the data, both at a hardware level 
as the data was being taken and in software afterwards. 
The hardware injections were performed to validate the 
end-to-end analysis, including parameter estimation on 
detection candidates, whereas the injections in software 
serve as a useful comparison, free of any calibration er- 
ror in the detectors. Here we report on the analysis using 
six hardware and software injections, including waveform 
models for binary neutron star (BNS), neutron star - 
black hole binary (NSBH) and binary black hole (BBH) 
simulations, described in section \LU\ One of these hard- 
ware injections was performed without the knowledge of 
the data analysis teams as part of the 'blind injection 
challenge'; it was successfully detected as reported in [7]. 
We use these injections to illustrate the possible implica- 



The Optimal Network SNR is defined as 



tions for GW astronomy in section IV and we conclude 
in section IVl 



II. ANALYSIS 

A. Bayesian Inference 

Each data segment containing an injected signal was 
analyzed using a Bayesian parameter estimation pipeline 
to calculate the probability density function (PDF) of 
the unknown parameters of the waveform model. We 
will call 9 the vector containing those parameters. The 
actual content and dimension of 9, i.e. the dimensionality 
of the parameter space, depend on the waveform model 
used for the analysis (see section II B ) . 

The posterior distribution of 9 given a model H is given 
by Bayes' theorem. 



pi9\{d},H) 



p{9\H)pi{d}\9,H) 
Pi{d}\H) 



(1) 



where p{9\H) is the prior distribution of 9, describing 
knowledge about the parameters within a model H be- 
fore the data is analyzed, and p({d}\9, H) is the likelihood 
function^ denoting the probability under model H of ob- 
taining the dataset {c?} for a given parameter set 9. The 
likelihood is a function of the noise-weighted residuals af- 
ter subtracting the model from the data, and is thus a 
direct measure of the goodness of fit of the model to the 
data. 



SNR: 
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where the sum is taken over each detector det, with Sdet 
the signal in that detector and Sdetif) its noise power 
spectral density (PSD). 

Our model for the likelihood function is based on the 
assumption that the noise is stationary and Gaussian, 
and uncorrelated in different frequency bins. Although 
we do not expect this assumption to be precisely true for 
real detector noise, limited investigations suggest that 
this is an acceptable approximation when the data is of 
good quality [34] . 

The denominator of eq. [T] P{{d}\H) = Z^, is the 
evidence for the model H. As it is a normalization con- 
stant, the evidence does not affect the estimation of the 
parameters for a particular model i?, but it does allow 
us to compare the ability of different models to describe 
the data. The Bayes factor between two models, which 
quantifies the relative support given to each model by the 
data, is given by 



^^p{9\H,)p{{d}\lH,)d9 

P{{d}\H2) J^p{X\H2)pi{d}\X, H2)dX 



Zi _ P{{d}\Hi 

Z2 



(3) 



where Hi and H2 represent the two different models and 
9, A their respective model parameter vectors. The pa- 
rameter spaces can be completely unrelated, even of dif- 
ferent dimensionality, which allows us to compare mod- 
els which range from 9 parameters in the case of non- 
spinning systems, up to 15 in the case of a system with 
two spinning components. Throughout we will quote ev- 
idence relative to the Gaussian noise model, ln(Z) = 

\n{ZH) - In(ZGaussian)- 

The high dimensionality of the parameter space and 
the complicated structure of the likelihood function make 
it impractical to exhaustively calculate posterior quanti- 
ties. Instead we rely on stochastic sampling of the poste- 
rior PDF, which provides us with an approximation of the 
underlying true distributions, e.g. by binning the sam- 
ples to produce histograms. The analysis was performed 
using LALInference [33], which allows calculation of the 
prior, likelihood and templates using standardised func- 
tions. Because accurate sampling of the posterior PDF 
is a difficult task, we used two independent sampling al- 
gorithms on all but a few cases, cross-comparing results 
to confirm convergence. These were based on Markov 
Chain Monte Carlo (MCMC) [33 |3S] and Nested Sam- 
pling [37j techniques, and were included in the LALIn- 
ference package. We found good agreement between the 
two techniques giving us confidence in the results. For 
a sample of the cases we used a third algorithm, Multi- 
Nest [38], as an additional check. The Nested Sampling 
and MultiNest algorithms directly produce an estimate 
of the evidence |39| ; we also computed the evidence from 
MCMC results using direct integration [?D]. 
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B. Waveforms and source parameters 

Unless otherwise specified we shall use a system of 
units in which c ~ G ^ 1. The waveform as measured at 
a generic detector can be written in the frequency domain 
as [H]: 

sU.e) = [F+{a,5,i^)h+{},e') (4) 

where: 

• F^{a, S, ip) and F-^ (a, (5, if)) are the known antenna 
beam pattern of the detector, which gives the am- 
plitude response of the antenna to the + and x 
polarizations. 

• a and 5 are the right-ascension and declination of 
the source; 

• -0 is the polarization angle (see, e.g., |42|): 

• At is the time delay between the arrival of the 
signal at the detector and at a common reference 
frame (e.g., the center of the Earth), a known func- 
tion of a and 5] 

• /i+(/, 9') and /ix (/, 9') are the two independent po- 
larizations of the signal, with 9' = 9 \ {a, 5, ip} 

The actual form of hj^{f, 9') and /ix (/, 9') depends on 
which model one is considering. 

The GW signal can be written using the post- 
Newtonian (pN) expansion [43] . We used a total of 
nine different waveform models (or approximants) in 
our analysis, including the TaylorF2 (which gives a 
pN-based, analytic expression for the frequency- domain 
waveform) [44l |45] . SpinTaylorT4 (which is a time- 
domain pN model supplemented with spin precession 
equations) [15] and IMRPhenomB (which is a TaylorF2- 
like waveform model supplemented with a merger- 
ringdown and calibrated to agree with numerical rela- 
tivity simulations) Tf approximants with binary spin 
parameters either disabled, set aligned to the orbital an- 
gular momentum, or enabled and allowed full freedom of 
orientation. Table |l] lists the approximants used. For ex- 
ample, the TaylorF2 waveform at Newtonian amplitude 
order can be written in the frequency domain combined 
with the stationary phase approximation j41| : 

h.(f) - Kf) (5) 

hAD^^'^Kf) (6) 



where: 

M/) = ^/Jr^/^e-^W) (7) 

with V = (7r(7Tii + m2)j) ^ . 

We see that the two polarizations depend on six pa- 
rameters (the ttfc are known functions of M. and r], [44] . 
up to iV = 7) : 

• M. The chirp mass of the system: M = 
[mim2)^^^{mi + m2)~^^^, where mi and m2 are 
the component masses; 

• rj. The symmetric mass ratio, defined as: t] — 

(mim2)/(mi + TO2)^; 

• t. The inclination angle, i.e. the angle between the 
orbital angular momentum vector and the line of 
sight; 

• c?i. The luminosity distance to the source; 

• tc- An arbitrary reference time, usually chosen to 
be the time of coalescence of the binary; 

• (jjc- The phase of the waveform at the reference 
time tj,. 

Together with a, S and tp, these six parameters form 
the 9-dimensional vector 9 for the TaylorF2 model. The 
upper value of the index k in eq. ]9] determines the phase 
pN order, with = being the lowest order (OpN) and 
N — 7 the 3.5pN order. We assume that the eccentricity 
of the binary system is negligible, as radiation reaction 
circularises orbits efficiently before radiation enters the 
detector's observational band [15]. We can also see from 
equations [5] and [6] that the inclination l and distance c?l 
are strongly correlated. This correlation is illustrated in 
section IIIII 

When spins are present and we constrain the spin vec- 
tors to be aligned with the orbital angular momentum 
of the binary, we need two additional parameters bring- 
ing to 11 parameters. We use the spin magnitudes ai 
and 02, defined as Ui = {si ■ L)/mf where si and rrii 
are the spin and mass of the object i, and L is the unit 
vector along the orbital angular momentum. To account 
for both aligned and anti-aligned case, we allow ai to 
be in range [—1, 1]. When generic spins are considered, 
we need six additional parameters, and 9 becomes a 15- 
dimensional vector. In this case the spin magnitudes a;, 
defined as = jsil/mf and in the range [0, 1], and four 
angles specify the orientations. 

The presence of arbitrary spins will alter both the am- 
plitude and phase of the signal. The pN expansion of the 
phase evolution is changed to include additional terms 
dependent on the components of the spin vectors. If the 
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spins are not aligned with the orbital angular momentum, 
then the spin-orbit and spin-spin coupling will cause the 
orbital angular momentum vector to precess around the 
total angular momentum vector. As this implies that t 
(and tp) become time dependent, we will use the value 
when the GW frequency is 40 Hz. We refer the reader 
to [35] for more details about SpinTaylorT4. 

The pN approximation is valid for the inspiral phase of 
the binary, when the two objects are still far apart [49l 
[50] . For this reason, the TaylorF2 and SpinTaylorT4 
waveforms must terminate at an approximate point 
where they begin to break down, with the TaylorF2 end- 
ing at the innermost stable circular orbit [HI [45] fre- 
quency of a test particle orbiting a Schwarzschild black 
hole, and the SpinTaylorT4 terminating at the minimum 
energy circular orbit [46 . IMRPhenomB waveforms ex- 
tend the signal into the merger and ring-down stages of 
the binary coalescence, making use of a phenomenologi- 
cal model tuned by comparison with numerical relativity 
simulations. These later stages take place at higher fre- 
quencies than the inspiral, hence are more important for 
higher-mass systems where the frequency scale is lower, 
and where the merger and ring-down contribute a larger 
proportion of SNR by being in a more sensitive band of 
the detectors. 

For the analysis of non-blind hardware and software 
injections (the first seven models of Table we con- 
sidered the pN expansion either up to the 2pN order 
in phase or up to the 3.5pN order in phase (2.5pN for 
spin effects when included), and Newtonian order (OpN) 
in amplitude. For the blind hardware injection we used 
approximants at 2.5pN order in phase, as the blind in- 
jectio n only included terms up to this order [7] [T^ (see 
Sect. IIIC below). The different waveforms use differ- 
ent post-Newtonian expansions of the binary phase, and 
do not give identical results in general. When analyz- 
ing a real GW signal, the ignorance of the higher order 
pN terms will cause a systematic bias in the inference of 
the binary parameters, and using a range of different ap- 
proximants here will give us an indication of the size of 
this error compared to the statistical uncertainty in the 
results. 



Model 


Name 


Spin 


Merger and 


Ref. 






effects 


ring-down 




TF2 


TaylorF2 


no 


no 


|44| 


TF2_2 


TaylorF2 Q 2pN 


no 


no 


|44| 


TF2_RS 


TaylorF2 RedSpin 


aligned 


no 


|45| 


ST_NS 


Non-spinning STPN 


no 


no 


|46| 


ST_SA 


Ahgned spin STPN 


aligned 


no 


|46| 


ST 


Full STPN 


yes 


no 


|46| 


IMRPB 


IMRPhenomB 


aligned 


yes 


|47| 


TF2_25F] 


TaylorF2 Q 2.5pN 


no 


no 


|44| 


ST_25p[ 


Full STPN @ 2.5pN 


yes 


no 


|46| 


Used only for the blind hardware injection 



TABLE I. List of the waveform models used for the analy- 
sis. 'Aligned' refers to both spin vectors being aligned to the 
orbital angular momentum of the binary. 'STPN' refers to 
the SpinTaylor post-Newtonian waveform for precessing bi- 
naries in LAL [33]. In all models, the inspiral phase of the 
binary evolution is described by the post-Newtonian (pN) ex- 
pansion to 3.5pN order (2.5pN for spin effects when included) 
in phase, and Newtonian order (OpN) in amplitude, unless 
otherwise specified. 



analyses using waveform models that account for possi- 
ble spins, the prior on the spin magnitudes, ai,a2, was 
taken to be uniform in the range [0, 1] (range [—1, 1] in 
the spin-aligned cases), and the spin angular momentum 
vectors were taken to be isotropic. 

The computational cost of the parameter estimation 
pipeline precludes us from running it on all times; there- 
fore, the parameter estimation analysis relies on an esti- 
mate of the coalescence time as provided by the detec- 
tion pipeline [7|. In practice, a 200 ms window centered 
on this trigger time is sufficient to guard against the un- 
certainty in the coalescence time estimates from the de- 
tection pipeline, see for instance [ST] [55] . Our results are 
not significantly affected by other astrophysically sensible 
choices of priors. For the SNRs used in this paper, our 
posteriors are much narrower than our priors for most 
parameters. 



Priors 



D. Data description 



As shown in eq. ^ the posterior distribution of 9 de- 
pends both on the likelihood and prior distributions of 
9. We used the same prior density functions (and range) 
for all the analyses, uniform in the component masses 
with the range IMq < mi, 2 < 30 M©, and with the to- 
tal mass constrained by mi -|- TO2 < 35 . This range 
encompasses the low-mass search range used in [7] , where 
1 M© < TOi,2 < 24 M0 and mi + 7712 < 25 M©. The prior 
density function on the location of the source was taken 
to be uniform in volume, constrained between luminosity 
distances € [1,100] Mpc. We used an isotropic prior 
on the orientation of the orbital plane of the binary. For 



Data from the multi-detector network consisting of two 
LIGO instruments (HI and LI) and Virgo (VI) were used 
coherently in the analysis. From each detector a total of 
32 seconds of data [-30, +2] s around the GPS time of 
the injection were analyzed. The initial frequency of the 
analysis was 40 Hz, which is low enough so that the de- 
tectors are not significantly sensitive to the start of the 
signals and template. The sampling rate was 2048 Hz 
(4096Hz for the BNS, section III A2), corresponding to a 
1024 Hz (2048 Hz) Nyquist frequency, high enough to in- 
clude the entire waveform for all models except IMRPhe- 
nomB, which has a negligible contribution to the signal- 
to-noise ratio at higher frequencies. The 32-second seg- 
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HW/ 


M 


mi 


7712 


dL 




l"i 


"2 


SNR 




SW 


(Me) 


(Mo) 


(Mo) 


(Mpc) 


n 








IIIAl 


HW 


3.865 


4.91 


4.02 


36.2 


26 








13 


III A 2 


HW 


1.502 


1.808 


1.647 


5.94 


138 








36 


IIIBl 


SW 


4.76 


6 


5 


30.0 


1.1 


0.6 


0.8 


19 


niB2 


SW 


2.99 


10.0 


1.4 


16.0 


0.5 


0.7 





13 


IIIC 


HW 


4.96 


24.81 


1.74 


24.37 


109 


0.57 


0.16 


16 



TABLE II. Parameters of hardware (HW) and software (SW) 
injected signals discussed in section |III[ Non-spinning injec- 
tions were generated using the EOBNR [53] waveform model, 
whereas spinning injections used the SpinTaylor model. The 
SNR column shows optimal network SNR, eq. [2] 





TF2 


TF2 


TF2 


TF2 


ST 


ST 


ST 


ST 


IMR 






_2 


_25 


_RS 


_NS 


_SA 




_25 


PB 


IIIAl 


69 


70 




70 


69 


71 


71 




73 


III A 2 


686 


699 




694 


685 


697 


694 




668 


IIIBl 


154 


146 




153 


158 


158 


157 




155 


IIIB2 


48 


48 




50 


50 


52 


64 




52 


IIIC 






136 










213 





TABLE III. Logarithm of the evidence \n{Z) relative to the 
Gaussian noise model for each injection and each waveform 
family obtained via direct integration. The numbers come 
with statistical error bars of ±5. 



ments of time-domain data were Tukey windowed, with 
a 0.4 s roll-off on either side of the segment. The PSD of 
the instrumental noise was estimated using 1024 s of data 
after the end of the analyzed segment. We verified that 
varying the methods of PSD estimation (using ±512 sec- 
onds spanning the signal trigger time, and using median 
and mean estimation methods), had a negligible effect on 
the parameter estimation results, smaller than the sys- 
tematic uncertainties of parameter estimation, although 
varying the time during which the PSD is estimated can 
have a significant effect (see section HID). 
Calibration 



errors, which can influence the recon- 
structed amplitude, phase, and timing of the data [50] . 
have the potential to affect parameter estimation results. 
An analysis of the effect of calibration errors, which con- 
sidered mock errors similar to those expected during the 
86 and VSR2/3 runs, concluded that such errors are un- 
likely to cause a significant deterioration in parameter 
estimation accuracy [32] at the moderate SNRs consid- 
ered in this paper. 



III. SIMULATIONS 

Over the course of the LIGO S6 and Virgo VSR2/3 
science runs a series of hardware injections have been 
carried out where the arm lengths of the three detectors 
were physically changed to simulate the passing of a GW. 

As an end-to-end test of the search pipeline during 



the science runs, a signal was added to the data via a 
hardware injection, without the data analysts' knowledge 
('blind'). The parameters and template family were re- 
vealed only after the search was complete, and we include 
a retrospective analysis of this blind injection in section 
III C (see [12] for parameter estimation carried out before 
the injection was revealed). 

We have also added several simulated software injec- 
tions into real detector noise from those runs. Below, we 
describe the results of parameter estimation analysis on 
several of these injections, whose parameters are listed 
in Table [llj We present the posterior probability den- 
sity functions on the source parameters using a coherent, 
multi-detector data model with the seven waveform fam- 
ilies described above. In Appendix |A] we present 90% 
credible intervals obtained using the waveform models in 
Table |T] for each injection in Table [III 



A. Non-spinning hardware injections 

1. Binary black hole 

We first describe the analysis of a hardware injection 
corresponding to a binary black hole with non-spinning 
4.91 M0 and 4.02 components, injected at a network 
SNR of 13 (SNRs of 8.7 in the Hanford and Livingston 
detectors and 4.4 in the Virgo detector). This injection 
was made with effective-one-body - numerical-relativity 
(EOBNR) waveforms, using the EOBNRvl version as 
described in [53| . 

Figure [T] shows a comparison of the posterior PDF of 
the mass parameters for different models. The chirp mass 
has a low statistical uncertainty of ^ 1 %, with the great- 
est statistical uncertainty claimed by models that allow 
for non-zero spin magnitudes due to inter-parameter de- 
generacies. The 90% credible intervals on the chirp mass 
obtained with TaylorF2 templates just exclude the true 
value, an indication of the systematic bias due to the 
waveform differences between the injected EOBNR wave- 
form and these templates. These differences are also re- 
sponsible for the particularly strong bias in the mass ratio 
for 2.0pN TaylorF2 templates, exceeding the statistical 
measurement uncertainty in this parameter. In contrast, 
we find that by using the 3.5pN TaylorF2 templates the 
systematic bias is reduced to less than the statistical un- 
certainty at the 90 % credible interval. The uncertainty 
in the mass ratio is typically larger than that in the chirp 
mass, leading to the characteristic thin, correlated joint 
distribution for mi and m2 evident in figure [l] (right). 

Figure [2] shows a selection of PDFs for the extrinsic 
parameters of the BBH source: the recovered sky posi- 
tion, distance and inclination angle. The binary can be 
constrained to two neighboring regions of the sky repre- 
senting the reflection of the source location through the 
plane of the detectors, which cuts this contiguous region 
in two. The correlations between intrinsic parameters 
(masses and spins, see |54H56j ) and extrinsic parameters 
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FIG. 1 . (left) Posterior probability distributions for the chirp mass M of the non-spinning BBH hardware injection (section 
III All for the seven signal models considered. The injected value is marked with a vertical red line, (right) Overlay of 90% 
probability regions for the joint posterior distribution on the component masses mi , m2 of the binary. The true value is marked 
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FIG. 2. Joint posterior probability regions for the location and inclination angle of the non-spinning BBH hardware injection 
(section |lll A 1 1 for the seven signal models considered, (left) The binary is constrained to two neighboring regions of the sky. 
(right) The distance and inclination, like the sky location, are estimated with a similar accuracy in models that include or 
exclude spins. 
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FIG. 3. Posterior probability distributions for the dimensionless spin magnitude of the heavier (left) and lighter (right) 
components of the binary from the non-spinning BBH hardware injection (section III A 1 1, as inferred in the model ST (table|T]l, 
full-spin STPN. The injection was made with — a2 = 0. 
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are generally weaker than correlations between different 
parameters within the same group. Thus both models 
that include and models that exclude non-zero spins have 
similar statistical measurement uncertainties for extrinsic 
parameters. 

Finally, figure |3] shows individual posteriors for the two 
dimensionless spin magnitudes, obtained using the model 
ST (table|l]), fuh -spin STPN waveforms. Although spin is 
not very strongly constrained, particularly for the lower- 
mass secondary, both spin measurements are consistent 
with the true value of spin. The absence of strong 
constraint is due to both the small difference in masses 
and nearly face-on inclination. This inability to constrain 
the spin is reflected in the evidence for each template 



family, shown in table III where all models have the same 
evidence within the error bars. 



2. Binary neutron star 



We also analyzed a hardware injection simulating a 
non-spinning neutron star binary system with 1.81 
and 1.65 M0 components, injected at a network SNR of 
36 (SNR of 26 in the Hanford detector, 17 in the Liv- 
ingston detector and 20 in the Virgo detector). This in- 
jection was made with the same EOBNR waveform fam- 
ily as in section [ill A 1| above. 

Figure [4] shows a comparison of the posterior PDF 
of the mass parameters for the different models. The 
chirp mass has a low statistical uncertainty of < 0.1% 
thanks to the high network SNR, with the greatest sta- 
tistical uncertainty claimed by models that allow for non- 
zero spin magnitudes due to inter-parameter degenera- 
cies, similar to section [III A 1| The greater statistical 
uncertainty of the IMRPhenomB waveform model stems 
from this model not describing such a low-mass system 
{nil + 1712 < 10 M0) very well |47j . C orrespondingly it 
has the lowest Bayes factor in table III The 90% credible 
intervals on the chirp mass obtained with TaylorF2 tem- 
plates include the true value. Unlike in section [III A 1 



the merger present in the injected EOBNR waveform 
happens at a frequency ( 1.2kHz) where the detectors 
are not very sensitive. Figure [5] shows the extrinsic pa- 
rameters of the BNS source: sky position, distance and 
inclination angle. In this case, the data were such that 
the injected distance and inclination angle lay far outside 
the 90 % credible intervals. This could be due to a par- 
ticularly unlikely noise realization, some non-stationarity 
or non-Gaussianity in the noise (i.e., a "glitch"), or some 
issue with the way the hardware injection was carried 
out. 



B. Spinning software injections 

1. Binary black hole 

We simulated the signal from a binary black hole with 
misaligned spinning components with a SpinTaylor soft- 
ware injection. This binary consists of 6M0 and SMq 
black holes, with dimensionless spin magnitudes of 0.6 
and 0.8, misaligned with the orbital angular momentum 
by angles of 40° and 150°, respectively, and an optimal 
network SNR of 19 (SNR of 17 in the Hanford detector, 
6.5 in the Livingston detector and 4.9 in the Virgo detec- 
tor). The misalignment between the spins and the orbital 
angular momentum causes the plane of the binary to pre- 
cess, producing both amplitude and phase modulations 
in the received GW signal (see section II B I . 

Figure [6] shows a comparison of the posterior PDFs 
of the mass parameters inferred with different template 
models. The mass ratio is again severely biased for the 
2pN TaylorF2 model, which is not surprising for a 3.5 pN 
injection with spinning components as both the pN order 
and the spin alter the phase evolution of the signal. 

Figure [7] shows the recovered sky location, distance and 
inclination angle of the source. For all those parameters 
the injection value lies just outside of the 90 % credible 
regions. 

The spin magnitudes are again poorly constrained, 
with non-zero support across the entire allowed range, 
as shown in figure [8j Even though the injected signal 
was simulated from a system with high spin magnitudes 
(0.6 and 0.8), the low difference in masses, the near anti- 
alignment of the spins with the orbital angular momen- 
tum and especially the face-on inclination conspire to give 
the poor spin estimates. In addition, with this weak pre- 
cession effect, the spin tilts (angles between the spin vec- 
tor and the orbital angular momentum) are also poorly 
constrained. A comparison of the evidences in table |III| 
indicates that all models have the same evidence within 
the error bars. 



2. Neutron star - black hole 

We simulated the signal from a neutron star - black 
hole system with misaligned black hole spin components 
with a SpinTaylor software injection. This binary con- 
sists of a 10 M0 black hole with dimensionless spin mag- 
nitude of 0.7, misaligned with the orbital angular mo- 
mentum by an angle of 130°, and a 1.4 Mq non-spinning 
neutron star, at an optimal network SNR of 13 (SNR of 
7.9 in the Hanford detector, 9.2 in the Livingston detec- 
tor and 3.6 in the Virgo detector). 

Figure |9] shows a comparison of the posterior PDFs 
of the mass parameters inferred with different template 
models. Figure [lO] shows the recovered sky location, 
distance and inclination angle of the source. The joint 
distance-inclination posterior degeneracy prevents either 
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FIG. 4. (left) Posterior probability distributions for the chirp mass A4 of the non-spinning BNS hardware injection (section 
|IIIA2[ ) for the seven signal models considered. The injected value is marked with a vertical red line, (right) Overlay of 90% 
probability regions for the joint posterior distribution on the component masses mi, m2 of the binary. The true value is marked 
by the blue star. 
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FIG. 5. Joint posterior probability regions for the location and inclination angle of the non-spinning BNS hardware injection 
(section III A 2 1 for the seven signal models considered, (left) The binary can be constrained to two regions of the sky 



representing the reflection of the source location through the plane of the detectors, (right) The distance and inclination, like 
the sky location, are estimated with a similar accuracy in models that include or exclude spins. The data were such that the 
injected distance and inclination angle are far outside the 90 % credible intervals. 



parameter from being measured precisely. The model de- 
scribing most accurately the data, with the highest Bayes 



factor in table III Full STPN, delivers the most accurate 



parameter estimates, as expected since this model was 
used for the injection. 



Figure [TT] shows the posterior PDFs for the spin pa- 
rameters. The spin magnitude of the neutron star is un- 
constrained due to the large difference in masses. The 
spin of the more massive black hole encodes more infor- 
mation in the waveform, but is also poorly constrained 
due to the almost face-on inclination. 



C. Blind hardware injection 

The search pipeline described in [7] identified a GW 
candidate occurring on 16 September 2010 at 06:42:23 
UTC. A Bayesian analysis was performed using the algo- 
rithms and implementations described above, where pa- 
rameter estimates varied significantly depending on the 
exact model used for the gravitational waveform. 

Following the completion of the analysis, the event was 
revealed to be a blind injection. Further investigation re- 
vealed several problems with the pre-un-blinding param- 
eter estimation: 

• The template signal included phase corrections 
only up to 2.5pN order, which is an outlier in the 
post-Newtonian expansion, see [S7]. At that time 
(before the blind injection was revealed as such) pa- 
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FIG. 6. (left) Posterior probability distributions for the chirp mass M of the spinning BBH software injection (section 
for the seven signal models considered. The injected value is marked with a vertical red line, (right) Overlay of 90% probability 
regions for the joint posterior distribution on the component masses mi, m2 of the binary. 



-40^^ 



-50' 



-60° T 



1 1 

* 








TaylorF2 @3.5PN 

Aligned spin STPN 

Full STPN 

TaylorF2 RedSpin 

IMRPhenomB 

TaylorF2 @2PN 
Non-spinning STPN 





















21h 



20h 



19h 

a 



18h 




180° 



FIG. 7. Joint posterior probability regions for the location and inclination angle of the spinning BBH software injection (section 
III B 1 1, (left) The binary's true location lies just outside of the 90% credible interval, (right) The degeneracy in distance and 
inchnation prevents either parameter from being accurately constrained individually. 
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FIG. 8. Posterior probability distributions for the dimensionless spin magnitude of the heavier (left) and lighter (right) 
components of the binary from the spinning BBH software injection (section |III B as inferred in the model ST (table |l|, 
full-spin STPN; the true values are shown with vertical red lines. 
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FIG. 9. (left) Posterior probability distributions for the chirp mass M of the spinning NSBH software injection (section 
for the seven signal models considered. The injected value is marked with a vertical red line, (right) Overlay of 90% pro 
regions for the joint posterior distribution on the component masses mi, m2 of the binary. 
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FIG. 10. Joint posterior probability regions for the location and inclination angle of the spinning NSBH software injection 
(section III B 2 1, (left) The binary is localized well on the sky. (right) In this case, the true value hes outside of the 90% credible 
interval of the joint distance-inclination marginalized probability density function. 
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FIG. 11. Posterior probability distributions for the dimensionless spin magnitude of the heavier (left) and lighter (right) 
components of the binary from the spinning NSBH software injection (section |lll B 2[ ), as inferred in the model ST (table |l|, 
full-spin STPN; the true values are shown with vertical red lines. 
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FIG. 12. (left) Posterior probability distributions for the chirp mass Ai of the blind injection (section III C I for signal models 
at 2.5pN. The injected value is marked with a vertical line, (right) Overlay of 90% probability regions for the joint posterior 
distribution on the component masses mi , m2 of the binary. The bias introduced by an analysis with a model which disallows 
spin is clear. 
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FIG. 13. Joint posterior probability regions for the location and inclination angle of the blind injection (section |III C \. (left) 
The sky location is constrained to several distinct regions lying along a half circle on the sky. (right) Again, the characteristic 
V-shape degeneracy in distance and inclination is evident. 
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FIG. 14. Posterior probability distributions for the dimensionless spin magnitude (left) and tilt angle (angle between the spin 
vector and the orbital angular momentum, measured at 40Hz, right) of the heavier component of the binary from the blind 
injection (section [ill C I, as inferred with model full-spin STPN at 2.5 pN order in phase (table |T]|; the true values are shown 
with vertical red lines. The large difference in masses allows for the spin of the massive component to be measured. On the 
other hand, the spin of the light component is unconstrained. 
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rameter estimation was not carried out with tem- 
plates at this order, leading to a significant bias in 
the mass ratio, and hence the component masses. 

• The signal to be injected in the Hanford and Liv- 
ingston sites had the wrong sign, making the signal 
incoherent between the LIGO and Virgo detectors. 
This caused strong biases in the estimated values 
of the extrinsic parameters. 

• Lastly, an injection software bug artificially set to 
zero one of the phase terms in the injected wave- 
form. 

We present here an analysis using templates at 2.5pN 
order in phase, performed after the un-blinding. We arti- 
ficially introduce a sign flip in the templates for Hanford 
and Livingston detectors in order to match the injected 
waveform. Figure [T2| shows a comparison of the posterior 
PDFs of the mass parameters for TaylorF2 and SpinTay- 
lor waveform models. The fully spinning, precessing anal- 
ysis with SpinTaylor templates at 2.5 pN order in phase 
(i.e., the same template as used in the injection, up to 
the injection software bug) centers on the correct masses 
for this neutron star - black hole system, rrii = 24.81 M0 
and m2 = 1.74 Mq. However, the systematic bias due to 
fixing spin to zero in the analysis of this spinning injec- 
tion is very significant, and leads to the wrong conclusion 
that a BBH system is observed if the TaylorF2 model is 
used. 

Figure [13] shows the recovered sky position, distance 
and inclination angle of the source. The sky location 
recovered is contained in several distinct regions, spread 
in an arc on the sky. This behavior is consistent with 
two detectors contributing the majority of the SNR 15 for 
this source (SNR of 11 in the Hanford detector, 9.8 in the 
Livingston detector and 4.1 in the Virgo detector). On 
top of this first-order constraint, spin projection effects 
can break the degeneracy, as discussed in [55]. 

The large difference in masses, mi = 24.81 M©, TO2 — 
1.74 M0, in combination with the high inclination allows 
for the magnitude of the spin of the massive component 
to be measured with an accuracy ^ 10%, as shown in 
figure [M] Meanwhile, the spin of the light component 
is unconstrained, and the posterior on its magnitude is 
consistent with the prior. Correspondingly, the evidences 
in table |III| favor the spinning model. 

A deeper analysis of the data around the blind hard- 
ware injection revealed some non-Gaussian behavior of 
the Livingston detector before and during the injected 
signal. From figure [15] it can be seen that the analysis 
of the Hanford detector data or Livingston detector data 
alone produces mass estimates consistent with the full 
network analysis using data from all three detectors, but 
the Livingston detector data cannot constrain the pa- 
rameters as precisely as Hanford detector data, showing 
multiple modes. This leads us to conclude that detec- 
tor glitches have the potential to reduce parameter esti- 
mation accuracy. However, the coherent analysis of the 
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FIG. 15. Overlay of 90% probability regions for the joint 
posterior distribution on the component masses mi, 7712 of 
the binary for the blind injection (section IIICl, as inferred 
with model ST_25 (table ]T]|, full-spin STPN(at 2.5 pN order 
in phase), using data from the Hanford detector only, the 
Livingston detector only or the whole Hanford-Livingston- 
Virgo network (HLV Network). 



data from the multi-detector network increases robust- 
ness against the effect of a glitch in a single instrument 
by requiring that the recovered parameters are consis- 
tent with all datasets. This effect could be mitigated 
by including such glitches in the model of the data, as 
suggested in [59] . 



D. The effects of different noise realizations 



To exemplify the effect of noise on the recovered PDFs, 



ten software injections (see figure 16 1 with parameters 
identica l to th e non-spinning BBH hardware injection 
(section III A 1 ) were created with injection times in suc- 
cession at 10-seconds intervals. With each injection at 
a different time, the estimated PSDs were also (slightly) 
different (see figure 17), refiecting the slow change of the 



PSD as a function of time. All the signals were both in- 
jected and recovered using the TaylorF2 signal model to 
eliminate systematic effects from using the wrong wave- 
form family. Figure [16] shows the joint distributions of 
the recovered component masses as well as the distance 
and inclination. The spread of the 90 % credible intervals 
for different injection times is comparable to the spread 
observed using different waveform models, as shown in 
figures [T] and [2] illustrating the dependence on the spe- 
cific model and realization of the noise used. Fixing the 
injection time for a series of different PSDs gives nar- 
rower PDFs than in figure [16] For the particular data 
segment considered here, the same is also true for a fixed 
PSD and a series of injection times but to a lesser extent. 
We plan to address this in the future by performing the 
analysis with a model in which the PSD is not fixed, but 
parametrized as a function of a suitable set of unknown 
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FIG. 16. Overlay of 90% probability regions for the joint posterior probability distribution of the component masses mi, m2 
of the binary (left) and the distance and inclination (right) for a series of signal injection times. Both sets of distributions show 
a spread comparable to the spread observed between different waveform models, as shown in figures [l] and |2] 
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FIG. 17. The plots show in red 10 overlaid noise noise power spectral densities (from left to right, LIGO-Hanford, LIGO- 
Livingston, and Virgo, respectively), each computed over 1024 sec of data, and separated by 10 sec. In black is the PSD used 
in Figs. [T] and [2] The PSDs were used in the analyses whose results are shown in Figs. |16|and |18| and provide an indication of 
the fluctuation of the noise in the instruments over a time-span of minutes. 




FIG. 18. (left) Posterior probability distributions for the chirp mass M of the non-spinning BBH injection for signal model 
TF2 (table |l|, TaylorF2 at 3.5 pN order, for the hardware injection (section III All and 3 software injections of the same 
parameters injected 100, 500 and 1000 seconds before the hardware injection (see section p[IID[ |. The injected value is marked 
with a vertical red line, (right) Overlay of 90% probability regions for the joint posterior distribution on the inclination and 
distance of the binary. Different realizations of the noise, rather than differences between hardware and software injections, are 
the likely reason for the variations in recovered PDFs. 
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parameters to be marginalised over. For the ten software 
injections the true injected parameter is found within the 
given 90 % credible interval in 90 % of the injections, as 
expected. Figure [T8| shows a direct comparison between a 
non-spinning BBH hardware injection and three software 
injections replicated with the same parameters. All are 
evaluated using the PSD of the hardware injection with 
the software injections at 100, 500 and 1000 seconds be- 
fore the hardware injection. Again, there are variations 
in the recovered PDFs but the hardware injection is not 
an outlier. The variations are caused by the inherent 
dependency of the analysis on the particular noise real- 
ization used. Statistical fluctuations appear to be able to 
dominate any systematic bias introduced by the method 
of performing the hardware injection, such as the actu- 
ation function used in modelling the impulse applied to 
the mirrors to produce the desired signal. 



E. Computational cost 

The choice of approximant has strong implications for 
the computational cost of the analysis. Frequency do- 
main templates such as TaylorF2 and IMRPhenomB gen- 
erally require much less computational time than the 
time domain SpinTaylorT4, as they do not require nu- 
merical solution of the orbital frequency and spin evolu- 
tion equations or a Fourier transform before computing 
the likelihood. The time domain SpinTaylorT4 can in- 
clude the full six spin parameters, allowing precession 
of the orbital plane, whereas the frequency domain tem- 
plates can only include spin aligned with the orbit. The 
computational cost is also strongly influenced by the 
length of the templates, which ranged from 2.16 s for the 
blind injection (section III C I to IBs for the ENS hard- 
ware injection (section _ 



III A 2 ) , using a 40 Hz starting fre- 



quency. An additional effect is the higher frequency at 
which a ENS will merge, forcing us to analyze the data 
at a higher sampling rate. A higher signal-to-noise ratio 
also increases the run-time with more iterations required 
to reach the more contrasted posterior distribution. This 
effect can be somewhat managed by increasing the num- 
ber of parallel tempering chains in the MCMC algorithm 
[551 or by using parallel nested sampling runs, re- 
ducing the overall run-time at the expense of using more 
compute units in parallel. Other promising cost-saving 
techniques include waveform interpolation to ease the 
burden of computing millions of templates [e.g.. I60j. 

In a future triggered analysis, where the parameter 
estimation is run as a follow-up to a search algorithm, 
some guidance on the chirp mass from the search pipeline 
could be used to choose the length of the data segment 
to be analyzed so as to minimize the cost while allowing 
enough time for the full waveform to be generated. In the 
analysis presented here, the run-time varied from hours 
for the low SNR, short segment (high mass) injections 
using frequency-domain templates, up to several weeks 
for the high SNR, long segments (low mass) using time- 



domain fully precessing templates. 



IV. IMPLICATIONS FOR 
GRAVITATIONAL- WAVE ASTRONOMY 

Although the 5 injections in Table [IT] are not statis- 
tically sufficient to cover the range of source possibili- 
ties, we do have samples representing the three expected 
classes, ENS, NSEH and EEH, at a range of SNR and 
some including spin. This allows us to highlight some 
of the general features of parameter estimates that are 
likely to be encountered with future analyses of signal 
candidates in advanced gravitational-wave detectors. 

We have found, in accordance with expectations, that 
the chirp mass parameter can be resolved with a very 
small statistical uncertainty of only a few percent or even 
lower, even when the SNR is relatively low, as in the case 



of the EEH hardware injection |III A 1| and the NSEH 
software injection |III E 2[ which had network SNRs of 13 
(see figures [T] and |9|. With the lower mass and higher 
SNR ENS hardware injection, the systematic differences 
between waveform models dominate the statistical uncer- 
tainty on the chirp mass, as the signal accumulates more 
cycles in the sensitive band of the detectors. 

Previous parameter estimation studies have mainly fo- 
cused on estimating the typical size of the errors, primar- 
ily using the Fisher information matrix |61l . The inverse 
of the Fisher matrix provides lower bounds for the covari- 
ance matrix of maximum likelihood estimators (Cramer- 
Rao lower bound) [29 . In [62 a single detector analysis 
of non-spinning signals suggested that the chirp mass M 
could be estimated by LIGO with a fraction of percent 
accuracy and the mass ratio r] with an accuracy of a few 
percent. We find in our analysis (section |III[ ) that for 
the 5 injections considered this accuracy is reached only 



for the high SNR example in section HI A 2 In general, 
the full analysis recovers chirp mass M and the mass 
ratio rj at the percent and tens of percent levels, respec- 
tively. In addition, [52] underlined the dependence of the 
accuracy on the phase order of the post-Newtonian ex- 
pansion, which we illustrate in the figures of section [Hi] 
Multiple detector studies, based on both Fisher matrix 
estimates and numerical simulations [631 164| showed the 
difference of accuracy in parameter estimation, between 
intrinsic parameters (e.g., masses) and extrinsic parame- 
ters (e.g., sky position, the strongly correlated inclination 
and distance). The estimation of the extrinsic parame- 
ters mainly depends on the number and positions of the 
GW detectors, while the accuracy for intrinsic parame- 
ters depends mainly on the SNR, and they are nearly 
equally well estimated with a single interferometer, for 
the same SNR. The Fisher matrix calculation has known 
limitations [65] , and tends to underestimate the errors at 
low SNR [55], which our results confirm. 

In all examples, the mass ratio parameter is strongly 
affected by systematic differences between the waveform 
models. For those models which included spin, the mass 
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ratio parameter had a larger uncertainty due to correla- 
tions between the spins and mass ratio (see [54H56] ). but 
these broader posteriors did encompass the true values. 
This indicates that even for sources which are expected 
to have an insignificant spin, i.e. BNS systems, a con- 
servative analysis with fully spinning templates is desir- 
able. For systems involving black holes, performing the 
full spinning parameter estimation is essential, and fur- 
ther improvements in waveform modelling, informed by 
numerical relativity simulations, are crucial for reducing 
systematic biases and therefore providing more precise 
mass and spin estimates. 

For the extrinsic parameters, of which the most in- 
teresting are the source location and distance, the accu- 
racy of inference is limited by the statistical uncertainty 
coming from the finite SNR of the source. The over- 
all signal amplitude in each detector may be measured 
with a fractional accuracy ~ 1/SNR, but the amplitude 
in each detector is a complicated function of the loca- 
tion, distance and polarization angle (see eqn. |4|, which 
produces the complicated credible interval shapes seen 
above. Most notably, the precision to which the distance 
can be measured is limited by the presence of correlation 
with the inclination angle, which is true for any detec- 
tor network. On the other hand, the degeneracy in sky 
position can be greatly reduced by observing the signal 
strongly in multiple detectors; in the case of the three- 
detector network used here the sky position is restricted 
to two locations, symmetric when reflected across the 
plane containing the detectors, and with opposite incli- 
nations (see fig. [5]). This degeneracy is broken in the 
case of the spinning BBH signal ( III B 1| where the in- 
clination angle can be determined from the precession of 
the orbital plane, eliminating one of the two sky loca- 
tions (see figure [7]). As the size of the two sky locations 
is primarily determined by the location of the detectors 
and the SNR, these results are expected to be qualita- 
tively similar for the advanced detector network. When 
the signal amplitude is relatively low in one detector (see 
figure 10), more of the ring- like structure is observable. 



and the location cannot be well-determined. 



cross-comparison among different sampling algorithms. 

Although there has been rapid progress in the field of 
GW parameter estimation, a number of key questions re- 
main. Some of these are already being addressed, or will 
be addressed over the next few years before advanced de- 
tectors come online. Some of the most pressing questions 
are: 

• How precisely will we be able to measure spin mag- 
nitudes and spin tilt angles in precessing NSBH and 
BBH systems? 

• How important are systematic waveform biases due 
to imperfect waveform knowledge for various sys- 
tem classes? What are the accuracy requirements 
on waveform families for parameter estimation? 

• What is the best way to handle the long-duration 
signals in advanced detectors? Their extended 
bandwidth can cause current implementations to 
increase in runtime by factors of twenty-five or 
more. 

• How accurately can we measure finite-size and 
tidal-dissipation effects for systems involving neu- 
tron stars? 

• How is parameter estimation affected when the 
background noise is not stationary and Gaussian? 

We anticipate that the analysis methods used in this 
paper will be directly applicable to the first detections 
from the advanced detector era, allowing us to carry out 
astrophysical inference and to measure masses, spins, and 
sky locations of coalescing compact binaries. This, in 
turn, will improve searches for electromagnetic counter- 
parts, and ultimately allow us to solve a key inverse prob- 
lem of GW astrophysics: to reconstruct binary evolution 
and dynamical binary formation from the observed distri- 
butions of masses and spins of merging compact binaries. 
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Appendix A: Credible Intervals 



In Tables [Wl [Vj |VT) |VlTl [Vml we report the median 
value and 90% credible intervals for the one-dimensional 
posteriors in Ai, ij, mi, m2, and (11 obtained in the con- 
text of the different waveform models for the injections 
listed in Table |llj Our reported credible intervals are 
symmetric, excluding an equal posterior mass (in this 
case 0.05) below and above the interval: a symmetric 
(100 — a)% credible interval [6'inin, 6'max] satisfies 
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TABLE VII. Median value and 90% credible intervals on 
chirp mass (M), symmetric mass ratio (rj), masses (mi and 
m2) and distance (d_t) obtained in the context of various 
waveform models (see Table |l| for the injection discussed in 
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TABLE IV. Median value and 90% credible intervals on chirp 
mass {A4), symmetric mass ratio (rj), masses (mi and m2) 
and distance (dL) obtained in the context of various waveform 
models (see Table II| for the injection discussed in § IIIAl 
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TABLE VIII. Median value and 90% credible intervals on 
chirp mass (M), symmetric mass ratio (??), masses (mi and 
m2) and distance (dz,) obtained in the context of various wave- 
form models (see Table[l| for the injection discussed in § III C 
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TABLE V. Median value and 90% credible intervals on chirp 
mass (A^), symmetric mass ratio (r;), masses (mi and m2) 
and distance (di) obtained in the context of various waveform 
models (see Table II| for the injection discussed in § III A 2 
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TABLE VI. Median value and 90% credible intervals on chirp 
mass (Af), symmetric mass ratio (77), masses (mi and m2) 
and distance (di) obtained in the context of various waveform 
models (see Table |l| for the injection discussed in § IIIB 1 
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